Coulomb phase diagnostics as function of temperature, interaction range and disorder 
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The emergent gauge field characteristic of the Coulomb phase of spin ice betrays its existence via 
pinch points in the spin structure factor S in reciprocal space which takes the form of a transverse 
projector V at low temperature: S{q) ^ V ^ q\/ . We develop a theory which establishes the fate 
of the pinch points at low and high temperature, for hard and soft spins, for short- and long-ranged 
(dipolar) interactions, as well as in the presence of disorder. We find that their detailed shape can 
be used to read off the relative sizes of entropic and magnetic Coulomb interactions of monopoles 
in spin ice, and we resolve the question why pinch points have been experimentally observed for 
Hoi.7Yo.3Ti207 even at high temperature in the presence of strong disorder. 

PACS numbers: 75.10.Jm 05.30.Jp 71.27.+a 



Motivation: The study of 'topological' states of mat- 
ter has been one of the central enterprises of condensed 
matter physics over the years, encompassing as they do 
fractional quantum Hall states, spin liquids or topological 
insulators. They are distinct from conventional phases 
exhibiting local order - e.g. a ferromagnet or charge- 
density wave - but also from simple disordered phases. 

The central experimental challenge is to diagnose their 
existence in more than a negative way, such as the ab- 
sence of order. Especially for states with mobile charged 
excitations, non-local probes such as conductivity^ or 
interferometry^ have been used to detect edge states 
or fractionalised excitations; in their absence, diagnos- 
tics have often resorted to general thermodynamic argu- 
ments, such as the presence of power laws in thermody- 
namic quantities^ which while highly suggestive have no 
direct link to the nature of the effective degrees of free- 
dom. 

In the case of spin ice, a three-dimensional frustrated 
magnetic material on the pyrochlore lattice consisting of 
corner-sharing tetrahedra which exhibits a topological 
Coulomb phase, a pinch-point motif in the spin struc- 
ture factor (Figllja)) has been identified as a diagnostic 
of the emergent gauge field characteristic of that phase, 
and analysed in considerable detail^^^: it emerges as a 
consequence of a local constraint, the "ice rule", imposed 
by the interactions, that each tetrahedron have vanish- 
ing total (pseudo-)spin = XlieKI ~ O^i^iii Excited 
tetrahedra violating this constraint carry an emergent 
gauge charge^i^ii^i — and interact through both an en- 
tropic and a magnetic (energetic) -^^i^^ Coulomb interac- 
tion. 

It turns out that the pinch points persist to regimes far 
away from the Coulomb phase, most strikingly to high 
temperatures^^ even in the presence of strong dilution 
of the magnetic Holmium ions by non-magnetic Yttrium 
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code much useful and interesting information. A com- 
bined analysis of the correlations along the two axes of 
the pinch point which we call nodal^ and antinodal (see 
Fig[T]) reveals this information, in particular the disconti- 
nuity where they cross at the pinch point itself, as well as 
the nature of the vanishing of the nodal correlators, the 
importance of which has been emphasized by Henle}*^^ 
and Bramwell.^^ We find the possibility of qualitatively 
distinct correlation functions for gauge charges and spin 
correlators: the latter are enhanced, while the former 
are screened, as long-ranged interactions are introduced. 
The qualitative enhancement is weakened, but does not 
disappear, at high temperature, even in the presence of 
disorder; while at low T, one can read off the relative size 
of entropic and magnetic charges. 

Our analysis thus contributes not only to the study of 
diagnostics of topological phases, as well as their inter- 
play of disorder, but also has connections to the unusual 
correlations induced by dipolar interactions, which are 
also of much current interest in the study of cold atoms 

The core of our analysis consists of a formulation of 
the large-n expression for correlations of the soft-spin 
model^ in terms of quantities natural to the emergent 
gauge degrees of freedom, which then straightforwardly 
generalises to the other settings. This may be of more 
general methodological interest. Here it provides an at- 
tractively simple yet complete picture in terms of (un- 
usual forms of) Debye screening. 

The analysis starts from the Hamiltonian^ 
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where ^ and D parametrise the strengths of nearest- 
ions This calls for a theory determining under which neighbour exchange and dipolar interactions, respec- 



conditions pinch points provide what information. 

We hence devise a theory for the structure factor across 
a range of settings (temperature, nature of spin, range 
of interaction, level of disorder) and find that they en- 



tively and Vij = Tij + 2(5^j, where F is the adjacency 
matrix of the pyrochlore lattice. We will also have occa- 
sion to consider an 'ideal' dipolar interaction, which dif- 
fers from Eq. [T]only by terms vanishing as 0(r~-^) with 
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FIG. 1. (a) <S(q) (see Eq[5)) calculated in the [hhk] plane for 
the ice states at T = (b) <S(q) calculated for Ho2-ccYccTi2 07 
at T = 10 K and x = 0.3 (c) Diffuse scattering along nodal 
direction (dotted line in (a) and (b)) and perpendicular to it 
(antinodal direction, bold line in (a) and (b)) for both the 
cases near the pinch point (0, 0,47r). q is measured relative 
to the pinch point. Note the weak pinch point discontinu- 
ity (indicated by arrows) even at high T and disorder, (d) 
Experimentally observed diffuse scattering from Ref [l3 for 
Ho2-xYa,Ti207 at T = 10 K and X = 0.3. In the three dif- 
fuse scattering plots, qx G [— 47r, 47r] and qz G [— 67r, Gtt] and 
measured in units of a~l- . 



distance, the matrix form of which corresponds to the 
projector V onto the ground states of the ice model. 

The pyrochlore lattice is made of up and down tetra- 
hedra labelled by rj = ±1, each up tetrahedron meeting 
down tetrahedra at its four vertices and vice-versa. It 
can be viewed as a fee Bravais lattice with an up tetra- 
hedron (four-site basis) at each fee site. The unit cell 
is then a conventional cubic cell of side acubic = 2>/2a, 
where a is the nearest neighbour distance. Si are (hard 
or soft) spin variables pointing along local easy axes in 
the a = 1 ... 4 [111] directions pointing from the centers 
of the up tetrahedra to the vertices. The centers of all 
tetrahedra form a bipartite diamond lattice with bond 
length ad = V3acubic/4. 

Nearest-neighbour soft spins: The large-n approach^^ 
works demonstrably well for the nearest-neighbor (NN) 
Heisenberg (n = 3) model at low T^^^ and for the 
Ising (n = 1) model at T = 0.^ The solution requires the 
eigenvalues and eigenvectors of Vij. In momentum space, 
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There are thus two fiat bands at zero energy and two 



bands (one gapless, quadratic acoustic and one optical) 
with energies 1^3, 4 (q) = ( J/3)'U3^4(q) for the NN model 
(I^ = 0). The eigenvectors in the flat bands satisfy the ice 
rules, i.e. the sum of spins on each tetrahedron vanishes. 

The spin correlations are expressed in the diagonal ba- 
sis of V^(q) using a unitary transformation U{c\). Then, 
(S'c(q)S'/3(-q)) equals: 

E A(T) ^ g ^^M)UlM) i,,(q)^A(T)T 
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where the q dependence of U is suppressed for brevity in 
the second line in Eq[2j The length constraint {S^) = l/n 
per component is imposed by the Lagrange multiplier 

A(r). 

The re- writing in the second line in Eq[2]will help us 
interpret results at finite T in terms of (gauge) charges 
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equals the charge-charge correlation 



(screening) length. This can be seen by calculating 
(S'^(0)5'^(r)). For r ^ a^ this equals: 
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Note that the charge-charge correlation function is 
isotropic in space unlike the spin-spin^ correlation func- 
tion. Also, when the screening length is large, i.e. 
<^cubic^ the charge density {S'^) approaches the 

"bare" value and A ^ f (1 + 0.224112^^). In the 
case of soft spins, = T/(J/3) from equipartition and 
«Lbic«'=8A(T)T/(J/3). 

From Eq[2l we can also calculate the structure factor 
^(q) = ^c.(q) ^/3(-q))- Its pinch points at 

T = are due to the singularities in the U (q) matrix.— 
Concentrating on small (Qx^Qx^Qz) near the pinch point 
at (0,0, 47r) for clarity, 5(q) can be written as: 
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where we omit the contribution of the gapped optical 
branch, which is fully analytic. 

At T = 0, only the fiat bands contribute, leading to 
pinch points. Moreover, when g^c = 0, the structure fac- 
tor is precisely zero near the pinch point. At finite T, 
there is non-zero diffuse scattering in this nodal direc- 
tion,^ proportional to Q^, the width of which is controlled 
by <^cubic^7 the inverse screening length. In the perpen- 
dicular antinodal direction, the scattering arises from the 
fiat band states and is present even at T = (FiglTJc)). 
The discontinuity in 5(q) is washed out at finite T and 
the structure factor is analytic at the pinch point. In real 
space, spin correlations now decay exponentially with 
a correlation length ^/acubic = \\/ at low T. 
Thus, the spin correlations and the charge correlation 
decay with the same correlation length in the NN case. 
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The charge correlations at low T can be understood 
as Debye-Hiickel (DH) screenin^^ of charges S^. Un- 
usually, their size is not "quantised" (because of the 
continuous nature of the soft spins) and their Coulomb 
interaction is entropic. The charges effectively live on 
the dual diamond lattice sites which are the centers of 
the tetrahedra of the pyrochlore lattice. The varying 
number of spin configurations consistent with a given 
charge configuration leads to the entropic Coulomb in- 
teraction JCad ■ QiQj/rij where /C = nT / (2^/37r) and 

Qi = rji{SM)i' The screening length equals ^cubic^^D?/ = 
{8V^7rJCp{Q^))/T within standard DH theory. For soft 
spins, each tetrahedron hosts a charge of mean size 
Q at low T: indeed, setting p{Q'^) = T/(J/3) gives 

Nearest-neighbour (hard) Ising spins: Ising spins can- 
not simply be described in terms of modes with eigenval- 
ues v^{q) because the unit length constraint couples the 
modes non- linearly. However, the rewritten form, Eq. [21 
does generalise straightforwardly in the sense of an effec- 
tive theory. The low-T properties of spin ice can be de- 
scribed in terms of a dilute set of mobile, Coulombically 
interacting charges on the dual diamond lattice. The cru- 
cial difference to soft spins is that now these charges are 
quantised, and hence gapped. The parameters and 
<^cubic^ are fixed by the properties of these interacting 
charges, while the Lagrange multiplier A(T) enforces the 
spin length constraint on average. Note that Vj^{ci) are 
the eigenvalues of the Laplacian operator Tij — zSij on 
the pyrochlore lattice, where z is the coordination num- 
ber, (with a further shift of 86 ij to make the eigenvalues 
positive-definite) which turns into in continuum. 

We first consider the NN Ising model. At low T, the 
(exponential) majority of defects have = ±2 so that 
= ±4 may safely be ignored. The fraction of defective 
tetrahedra p is itself exponentially suppressed and equals 
2exp(-2J/3T). Since (5^) = 4p, this fixes = 4p. 
The Lagrange multiplier A(T = 0) = 1/2 since n = 1 
for Ising spins. The weight in the flat bands decreases 
according to at finite T, which we take into account 
by A(T) = 1/2 + a{T)p. Specifying ^nd A(T) fixes 
the inverse screening length ^^^^^^-^^/^g — 16p + 32a(T)/)2 
for the NN Ising model as we require the pinch points to 
vanish at any T 7^ 0. Q^(T) can be determined numeri- 
cally from the unit-length constraint and equals 0.44822 
as T ^ 0. The leading term in <^cubic^e exactly repro- 
duces the entropic screening length calculated in Ref [23I 
using DH theory where the gapped charges {Q = ±2) in- 
teract with an entropic Coulomb interaction of strength 
JC = T/2V37r. 

The pinch points now acquire a finite width in q space, 
determined by acubic^e? of order \q\ ~ for Ising 
spin^^ (see Figs[2];a),(d)) instead of - V^T^ for soft 
spins. This width can be observed in the half-width at 
half maximum (HWHM) of the diffuse scattering along 
the pinch point's nodal direction. In spin ice, the Ising 
spins point along the easy-axes and thus, the experimen- 
tally relevant structure factor 5(q) has additional phase 



factors,^ and equals 
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Ca ■ (e|| X q) e/3 • (e|| x q) 
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where neutrons are polarized along e\\ = -^(1,-1,0) 
and the (spin- flip) scattering measured along the e|| x q 
direction where q lies in the [hhk] plane. We display 
5(q) in the discussions of diffuse scattering rather than 
5(q). Note that in three dimensions, the DH screening 
length p~^l^ scales differently from the average distance 
between the thermal defects, ~ p~^^^. As in the soft 
spin case, the NN Ising case also has both the spin-spin 
correlations and the charge-charge correlations decaying 
with the same correlation length ~ at finite (low) 

T. 

Dipolar spin ice: We now turn to dipolar spin ice 
where long-ranged interactions between the spins, D ^ 0^ 
play an important role. Due to the dipolar interactions, 
the charges ^ now have an additional magnetic 
Coulomb interactionJ^ This can be conveniently seen 
by using the dumbbell modeU^ where each spin (point 
dipole) is replaced by an extended dipole of size ad placed 
on each bond of the dual diamond lattice. Then, upto 
0(l/rf ) terms, the interactions can be re- written as 



3\/3 ^ 



A^(Q./2)2 
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where A = ^ + | (l + ^f) D. 

Thus, the charges interact with both an entropic and a 
magnetic Coulomb interaction and the total coupling is 
JC = T/{2VS7r)^2V2D/{3V3). The magnetic Coulombic 
DH screening length equals <^cubic^M ~ ^^^^ pD. Here, 
p needs to be calculated self-consistently^^ and equals 
2exp(— A/T) as T ^ 0. Thus, the fraction of defective 
tetrahedra now depends on both J and D. 

The parameters of the generalized equation are now 
Q2 = 4p, A(T) = 1/2 + a{T)p and = t^l^ where 
a(T) is determined from the unit length constraint. The 
parameters go smoothly to the NN Ising values when 
J ^ D. The charge-charge correlators calculated from 
this generalized equation have the correct DH form when 
both interactions are present. 

Screening is now enhanced as the magnetic Coulomb 
interaction does not vanish with T. The contribution 
of the acoustic band is then suppressed by 0{T/D) for 
wave- vectors \q\ <C tZM near the pinch point as T ^ 0, 
because p/t^\f ~ T/D. Thus, the pinch points sur- 
vive but are weakened. More specifically, the disconti- 
nuity at the pinch points, when approaching from the 
nodal and the antinodal directions respectively, scales as 
1 - (1 + 4V^7rD/3T)-i at low T. Secondly, the HWHM 
for the scattering in the nodal direction at the pinch point 



scales as y (1 + 4V^7rI^/3T)y^. Both the discontinuity 
at the pinch point and the HWHM in the nodal direction 
show that the gauge charges interact both energetically 
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FIG. 2. (a) Diffuse scattering <S(q) calculated along the nodal 
line near the pinch point (0,0, 47r) for both the NN Ising 
model and dipolar spin ice (using HTO parameters), q is 
measured relative to the pinch point. The dotted line shows 
the scattering for the NN Ising model using the value of p ob- 
tained at T = 0.65 K with HTO parameters. The dotted line 
meets the antinodal line (~ 2.66) at q = 0.(b) "Correlation 
length" ^ defined using the HWHM of the diffuse scattering 
calculated using HTO parameters (red points) in the nodal 
direction. Blue points shows the corresponding NN Ising re- 
sult using the p obtained with HTO parameters, (c) Diffuse 
scattering near the pinch point (0,0, Air) in the [hhk] plane 
for HTO parameters at T = 0.65 K. Note the discontinuity 
even at finite T. (d) Diffuse scattering calculated using the 
same p in the NN Ising case. The pinch point discontinuity is 
rounded off. Farther away from the pinch points, (c) and (d) 
agree. In both the scattering plots, qx/'^ir G [—0.05, +0.05] 
and q^/27v G [1.95,2.05]. 



and entropically: (4V^7rI^)/(3T), the ratio of the mag- 
netic and the entropic Coulomb couplings, can be probed 
directly in these quantities. 

Fig [2fa) shows the calculated diffuse scattering along 
the nodal direction at low T using Ho2Ti207 (HTO) pa- 
rameter^Si of D = 1.41 K and J = -1.56 K. FigEJc), 
shows the survival of the pinch point at finite T using 
HTO parameters at T = 0.65 K. A corresponding "cor- 
relation length" , which corresponds to that of the charge- 
charge correlation, can be extracted from the HWHM in 
the nodal direction (Fig[2fb)). In the same figure, we 
show the correlation length if the calculated defect den- 
sity p is used as the input parameter in the NN Ising case, 
where the charges only interact entropically. The growth 
of the correlation length is slower in the dipolar spin ice 
case because of the magnetic Coulomb screening. 

Thus, the picture that emerges is that the charge- 
charge correlations are again exponentially decaying 

in real space with a scale <^cubic*f~^ ~ ^1 + ^^^^ ^ P 
signaling the sparse defects of the Coulomb phase. 
Scattering along the nodal directions near the pinch 



points, driven entirely by these charges, reflects the 
screening of these charges which is now enhanced by 
the presence of the magnetic Coulomb interactions. 
However, the spin-spin correlations have a power-law 
tail which is dipolar and decays as 1/r^ at all tempera- 
tures^ signaling the presence of long-range interactions 
and not of a Coulomb phase. The last feature can 
already be seen from the soft-spin problem where the 
spins interact through the 'ideal' dipolar interaction^ 
Hij = Vij = Sij — {SiSj)T^o,n^i^ where {SiSj)T^o,n=i 
is calculated using Eq [2] at T = and n = 1. Then 
the spin-spin correlations (^Q;(q)*5'/3(— q)) become 

~ ^) (^a(q)^/3(-q))T=0,n=l + - 

with A(T) = riT-i+vi+n^T ^^ rj.^^ chargc-chargc correla- 
tions are however extremely short-ranged and non-zero 
only for two tetrahedra sharing a common site. 

This observation holds the key to explaining the pres- 
ence of well-defined pinch points at high temperature 
T = 10 K and strong chemical disorder in experiments on 
spin ice^*" diluted with non-magnetic Y, Ho2-a^Ya^Ti207, 
where the constraint = must be massively violated. 
These pinch points are entirely due to the form of the 
long-range nature of the dipolar interactions on the py- 
rochlore lattice, which mimic the spin-spin correlations 
at T = OJ^ This can be seen explicitly in the leading 
term of a high-temperature series expansion generalized 
to include disorder where the probability of a site being 
occupied by a magnetic Ho^+ ion is (1 — x/2): 



{S{n)S{r2)) 



{SirY) = (1 




(5(ri)5(r2))T=o,n=i 

5 - ^ \ _ \ 
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From this, we calculate 5(q) for Ho2-a^Ya^Ti207 at x = 
0.3 at a high temperature of T = 10 K (Fig[Hb)). The 
diffuse scattering pattern is similar to the experimen- 
tal one (Fig [Ifd)). Note that the discontinuity at the 
pinch point is weak in D/T^ and the "nodal" correla- 
tions are broad and large (Fig[Hc)). However, higher- 
order terms in 1 /T will not necessarily preserve the dom- 
inance of the pinch point features and instead lead to rich 
low-temperature physics including a dilution dependent 
"residual entropy" and the possibility of glassiness^^i^i 

Conclusions: Defects with gauge charge {S^ 7^ 0) in 
both NN and dipolar models can be generated by thermal 
excitations or chemical disorder. For D 7^ 0, these carry 
both entropic and magnetic Coulomb charges, while they 
carry only entropic charge in the NN model. This leads 
to a qualitative change in the spin correlations which 
can be observed in equilibruim experiments probing the 
pinch points. We have provided a simple and explicit 
computational and conceptual framework encompassing 
all these cases. 
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